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利用 测 地 观测 中 搭载 式 相 位 参考 模式 观测 
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摘要 : 基于 盖 亚 (Gaia) 天体 测 量 数据 建立 的 盖 亚 天 体 参 考 系 (Gaia Celestial Reference 
Frame, GCRF) 与 基于 其 长 基线 干涉 测量 ( Very Long Baseline Interferometry, VLBI) 建立 的 国际 
天 球 参 考 系 (International Celestial Reference Frame, ICRF) 具 有 一 定 的 系统 差异 ， 需 要 对 盖 亚 
的 数据 进行 外 部 的 独立 评估 ， 以 建立 光学 和 射电 参考 架 连接 。 国 际 上 对 参考 架 的 连接 普遍 采 
用 在 光学 和 射电 波段 都 有 辐射 的 天 体 ， 在 瞳 端 ( 星 等 C15) 0 3E 4E, EIER (C13) 
用 射电 星 。 由 于 普通 VLBI 观测 申请 批准 的 周期 长 ， 难 度 大 ， 因 此 ， 我 们 尝试 将 一 部 分 射电 
星 的 观测 转向 数量 较 多 的 测 地 观测 ， 在 测 地 观测 中 搭载 相位 参考 模式 来 观测 射电 星 。 为 了 测 
试 这 种 观测 的 可 行 性 ， 我 们 在 IVS(The International VLBI Service for Geodesy and Astrometry ) 
fe LBA ( Large Baseline Array) 的 测 地 观测 中 搭载 相位 参考 观测 ， 选 取 Del Lib 和 HR 1099 作为 
实验 的 目标 源 。 观 测 得 到 了 这 两 颗 射电 星 的 VLBI 图 像 、 流 量 和 单 历 元 坐标 ， 且 坐标 精度 高 
于 前 人 ; 同时 结合 前 人 的 数据 ， 对 这 两 颗 射 电 星 进行 了 视差 和 自行 拟 合 ， 结 果 与 前 人 和 盖 亚 
的 结果 相符 ， 且 部 分 精度 更 高 。 验 证 了 测 地 模式 中 搭载 相位 参考 观测 射电 星 的 可 行 性 ， 为 后 
续 丰 富 射 电 星 VLBI 天 体 测 量 样本 ， 连 接 光 学 与 射电 参考 架 黄 定 了 一 定 的 基础 。 
关键 词 : 甚 长 基线 干涉 测量 ; 射电 星 ; 天 体 测量 ; 自行 ; 视差 
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盖 亚 光学 卫星 最 新 数据 DR3(Data Release 3) CARA, 然而 ， 人 们 发 现 DR3 的 视差 随 所 处 天 
区 、 天 体 角 径 大 小 、 星 等 和 颜色 的 不 同 而 存在 一 定 的 系统 误差 ， 基 于 盖 亚 类 星体 的 数据 (理论 上 河 外 
源 视 差 应 为 0) 得 到 的 平均 视差 0 点 误差 约 为 -29 uas!3 ， 且 基于 盖 亚 天 体 测量 数据 建立 的 盖 亚 天 体 参 
考 系 与 基于 其 长 基线 干涉 测量 建立 的 国际 天 球 参考 系 有 微小 的 差异 ， 因 此 需要 对 羡 亚 的 数据 进行 
外 部 的 独立 评估 ， 以 期 建立 光学 和 射电 参考 架 连 接 。 

国际 上 对 参考 架 的 连接 普遍 采用 在 光学 和 射电 波段 都 有 辐射 的 天 体 ， 在 暗 端 ( 星 等 = 15 ) 使 用 类 
星体 ,在 亮 端 (G<13) 使 用 射电 星 。 尽 管 目前 可 用 于 参考 架 连 接 的 射电 星 ( 100 颗 左 右 ) 的 数目 远 小 于 
类 星体 ( 数 千 颗 ) ， 然 而 与 类 星体 相 比 ， 射 电 星 大 多 位 于 星 等 的 亮 端 ， 具 有 更 高 的 天 体 测 量 精度 ， 且 
不 存在 如 类 星体 大 尺度 喷 流 的 影响 汪 ” ， 因 而 更 适合 作为 光学 亮 端 和 射电 波段 天 球 参考 架 连 接 的 目 
标 。 目 前 可 用 于 参考 架 连 接 的 射电 星 中 约 40 颗 的 射电 天 体 测量 参数 (主要 包括 位 置 、 视 差 和 自行 ) 来 
自 于 上 址 纪 的 观测 ， 精 度 相 对 较 低 5 。 为 了 更 好 地 实现 光学 与 射电 波段 天 球 参考 架 的 连接 ， 必 须 在 
盖 亚 服役 期 间 增 加 具有 VLBI 天 体 测 量 参数 的 射电 星 数目 ， 从 而 避免 由 于 观测 历 元 不 一 致 引入 的 误 
差 。 因 此 ， 尺 快 开展 大 规模 的 射电 星 其 长 基线 干涉 测量 尤为 重要 。 
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由 于 普通 甚 长 基线 干涉 测量 观测 申请 批准 的 周期 长 ， 难 度 大 ， 大 批量 的 射电 星 观测 需要 耗费 较 长 
时 间 ， 因 此 可 以 将 一 部 分 射电 星 的 观测 转向 观测 频繁 、 数 量 较 多 的 大 地 测量 观测 ， 在 测 地 观测 的 空 际 
穿插 射电 星 的 相位 参考 观测 。 为 了 测试 测 地 模式 中 相位 参考 观测 对 于 流量 较 弱 的 射电 星 观 测 的 可 行 
性 ， 我 们 尝试 在 IVS 和 LBA 的 测 地 观测 中 加 入 一 些 扫描 来 观测 射电 星 ， 这 就 是 搭载 式 相位 参考 观测 。 
结合 测 地 观测 计划 ， 我 们 选取 在 观测 当天 各 台 站 均 可 见 ， 且 俯仰 较 合 理 的 射电 星 Del Lib 和 HR 1099 
作为 实验 的 目标 源 。 

Del Lib 全 称 8 Librae， 又 名 HR 5586, HD 132742 或 HIP 73473 ， 是 天 空中 第 3 Z£(m, 24.91) 的 
Algol HAE? 。 它 由 一 颗 温 度 较 高 的 B8 V 型 主 序 星 和 一 颗 温 度 较 低 且 质量 较 小 的 K2 IV 型 亚 巨 
星 组 成 ， 且 这 颗 较 小 的 亚 巨 星 充 满 了 它 一 侧 的 洛 希 办 ， 物 质 正 在 从 这 颗 亚 巨星 转移 到 主 序 星 '"。 该 
双星 系统 轨道 周期 约 为 2.33 天 "1 。 文 [8] 指 出 ， 该 Algol 系统 是 一 个 具有 长 周期 轨道 的 三 体系 统 ， 第 
3 个 成 分 为 Fl IV 型 星 。Del Lib 的 射电 光斑 爆发 似乎 存在 一 个 大 约 50 天 的 周期 "|。 

HR 1099 又 名 V711 Tauri, HD 22468 或 ADS 2644A, 文 [9] 和文 [10] 在 1976 年 各 自分 别 用 光谱 和 
光度 测量 观测 到 了 HR 1099， 并 且 确 定 其 为 RS CVn( RS Canum Venaticorum) 型 近 密 双星 ， 由 一 颗 KI 
IV 型 亚 巨 星 和 一 颗 G5 IV-V 型 主 序 星 组 成 ' 沾 ， 轨 道 周期 约 为 2. 84 R, EVI KI 型 星 比 另 一 颗 
更 活跃 ， 在 其 色 球 层 、 过 渡 区 和 日 园 中 具有 高 温 物 质 明亮 发 射 的 特征 '“'。 该 双星 系统 的 轨道 倾角 
大 约 为 30”'" ， 该 RS CVn 系统 也 是 一 个 具有 长 周期 轨道 的 三 体系 统 ， 第 3 个 成 分 ADS 2644B 角 距 离 
大 约 为 6"， 据 推测 对 RS CVn 系统 的 物理 影响 很 小 [2 。 


1 观测 安排 


我 们 的 观测 分 为 两 次 ， 分 别 为 IVS 和 LBA 的 测 地 观测 。 表 1 列 出 了 两 次 观测 的 代码 、 时 间 、 频 
率 设置 和 参与 台 站 等 信息 ， 表 2 列 出 了 两 次 观测 的 参考 源 信息 。 测 地 观测 一 般 为 了 减少 数据 量 而 仅 使 
用 右 旋 ， 所 以 第 3 节 中 数据 结果 图 像 也 仅 有 右 旋 偏振 。 


表 1 观测 配置 信息 
Table 1 Setup of the observations 


Project Date V ops Data rates Duration Stati d 
code ( yyyy-mm-dd) / GHz / Mbps /h m 
AOVO10 2016-07-27 8. 40 1 024 4 HO, IS, K1, KE, KM, KG, PA, SH, T6, 
TS, UR, WW, YG 
V515C 2018-07-19 8. 60 1 024 10 BD, HH, KM, SH, SV, ZC 


ik: 两 字母 代表 台 站 缩写 : Hartrao( HO), Ishioka( IS) , Kashim11(K1), Kunming( KM), Koganei( KG), Parkes (PA), 
Sheshan ( SH) , Tianma65(T6), Tsukub32(TS), Urumqi ( UR), Warkworthl2m( WW), Yarral2m( YG)), Badary( BD), 
Hartebeesthoek( HH) , Svetloe( SV) , Zelenchukskaya( ZC) 。 


R2 射电 星 和 它们 对 应 的 参考 源 ( 坐标 为 观测 纲要 文件 VEX 中 使 用 的 ) 


Table 2 Radio stars and their corresponding calibrators (coordinates used in VEX files) 


Star name Calibrator name R. A. (hms) Dec. (° ' ') Separation/ ( ° ) Project code 
Del Lib 1454-060 14 56 41. 39331 —06 17 43. 2601 2.47 AOVO10 
HR 1099 CTA 26 03 39 30. 9378 -01 46 35. 804 2. 46 V515C 


IVS 的 观测 代码 为 AOV010， 该 测 地 项 目 是 亚太 地 区 大 地 测量 及 天 体 测量 的 一 次 观测 ， 于 2015 年 
开始 !5] ， 主 要 参与 国 包括 澳大利亚 、 新 西 兰 、 日 本 和 中 国 。AOV010 的 观测 波段 在 S/X 波段 ， 为 了 
获得 更 高 的 空间 分 辨 率 ， 我 们 仅 使 用 X 波 段 的 8.4 GHz 数据 ， 单 通道 带宽 16 MHz， 其 中 台 站 Kashim11 
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和 Koganei 只 有 8 200 MHz 到 8 600 MHz 的 窄带 X 波段 ， 总 数据 记录 速率 为 1024 Mbps。 本 次 观测 的 主 
要 目的 除 使 用 测 地 模式 观测 一 组 天 体 测量 源 外 ， 还 搭载 了 对 于 目标 射电 星 Del Lib 和 其 参考 源 1454- 
060 的 相位 参考 模式 观测 ， 两 者 角 距 为 2. 47”， 来 测试 相位 参考 模式 能 否 观 测 到 较 弱 的 射电 星 。 我 们 
将 Del Lib 和 其 参考 源 1454-060 的 扫描 长 度 各 设置 为 2 min， 然 后 充分 利用 该 测 地 观测 的 时 间 安 排 空 
W, F Del Lib 的 相位 参考 扫描 分 为 3 部 分 ， 每 部 分 长 度 大 约 30 min， 两 部 分 间隔 至 少 为 1 h， 以 便 获 
得 更 好 的 UV A, R Del Lib 的 总 时 间 约 为 36 min， 时 间 跨 度 约 为 4h。 

LBA 的 观测 代码 为 VS15C， 也 是 测 地 类 观测 ， 搭 载 了 针对 射电 星 HR1099 的 相位 参考 观测 XU 
频率 为 SAX 波段 。 我 们 仅 使 用 X 波段 的 8.6 GHz 数据 ， 单 通道 带宽 16 MHz, 2 比特 量化 ， 总 数据 记 
录 速 率 为 1 024 Mbps。 我 们 使 用 强 射 电源 CTA 26 作为 参考 源 ， 目 标 源 HR 1099 和 其 参考 源 之 间 的 角 
FEX 2.46°, HR 1099 的 扫描 长 度 为 180s，CTA 26 的 扫描 长 度 为 60 s。 本 次 观测 虽然 是 LBA 组 织 的 ， 
但 参加 台 站 除了 中 国 的 云南 昆明 站 40 m 天线 (KM) 和 上 海 余 山 站 25 m 天 线 (SH) 外， 还 有 俄罗斯 应 用 
天 文学 研究 所 的 BD, SC 和 ZC 天 线 ， 以 及 南非 的 HartRAO 26 m 天 线 (HH) ， 基 线 长 度 远大 于 仅 在 澳 
大 利 亚 本 土 的 观测 ， 空 间 分 辩 率 相 较 于 常规 LBA 更 高 。HR 1099 的 时 间 ( On-source) 约 为 75 min， 总 
观测 时 间 也 分 为 3 部 分 ， 时 间 跨 度 约 为 10 h, 


2 数据 处 理 


TVS 的 观测 数据 由 上 海天 文 台 相关 处 理 中 心 59 进行 相关 处 理 ，LBA 数据 由 ATNF( Australia Telescope 
National Facility) 的 DiFX 相关 机 "中 进行 相关 处 理 。 首 先 我 们 将 相关 后 的 数据 装载 进 VLBI 数据 校准 
软件 AIPS( Astronomical Image Processing System! ) , 

一 般 情况 下 测 地 观测 同时 使 用 S/X 两 个 波段 ， 为 了 获得 更 高 的 空间 分 别 率 ， 我 们 仅 需 要 频率 较 
高 且 信 噪 比 较 高 的 X 波段 数据 ， 因 此 需要 运行 AIPS 脚本 vibautil, (EFA vlbafqs 命令 将 S 和 X 波段 的 
数据 分 制 ， 仅 保留 X 波 段 数 据 ; 然后 运行 ANTAB 任务 载 人 台 站 antab 文件 ， 进 行 幅度 校准 。 

经 过 星 位 角 改 正 、 电 离 层 改正 和 地 球 定向 参数 (Earth Orientation Parameters, EOP) 改正 后 ， 对 于 

AOV010， 我 们 对 该 次 观测 中 各 台 站 都 有 条 纹 、 流 量 较 强 的 一 颗 源 0454-234， 求 解 仪器 时 延 ; 然后 对 
相位 参考 观测 校准 源 1454-060 进行 全 局 条 纹 拟 合 ， 将 解 代 入 目标 源 Del Lib。 在 此 步骤 中 ， 我 们 发 现 
由 于 昆明 站 (KM) 相位 跳动 非常 大 (可 能 是 观测 时 台 站 设备 出 现 相 位 或 时 间 跳 动 ' ) ， 导 致 条 纹 拟 合 
求解 运算 时 失败 的 解 非常 多 ， 经 过 多 次 尝试 后 ， 我 们 将 此 步骤 分 为 两 步 : (1) 在 去 除 昆 明 站 的 情况 下 
对 其 他 所 有 台 站 进行 条 纹 拟 合 ; (2) 在 上 一 步 解 的 基础 上 ， 仅 对 昆明 站 进行 条 纹 拟 合 ， 然 后 将 两 步 的 
解 亚 加 应 用 到 目标 源 Del Lib 上 ， 最 后 导出 校准 完毕 的 数据 。 
与 AOV010 初始 处 理 步 骤 类 似 ， 在 完成 星 位 角 、 电 离 层 和 地 球 定向 参数 改正 后 ， 对 于 V515C XL 
测 ， 首 先 对 校准 源 CTA 26 进行 全 局 条 纹 拟 合 ， 将 解 应 用 到 目标 源 HR 1099， 然 后 将 AIPS 校准 完毕 后 
的 数据 ， 导 入 VLBI 图 像 处 理 软件 DIFMAPPU ， 对 校准 源 CTA 26 进行 成 像 。 图 像 显 示 它 并 非 点 源 ， 
具有 了 明显 单 侧 喷 流 结构 ， 见 图 1(a) 。 所 以 我 们 将 DIFMAP 生成 的 图 像 模型 FITS 文件 ， 导 入 AIPS, 使 
用 AIS 任务 CALIB 进行 两 次 迭代 ， 以 减少 校准 源 复杂 结构 的 影响 (Del Lib 的 校准 源 1454-060 的 图 
像 显 示 其 结构 较为 紧凑 ， 见 图 1(b) ， 不 需要 此 步骤 ) ，AIPS 校准 完成 ， 导 出 文件 。 

如 果 台 站 在 观测 时 能 够 正常 测量 系统 温度 和 增益 ， 则 观测 完毕 后 ， 台 站 会 生成 各 自 的 antabfs X 
件 ， 然 后 上 传 至 相关 处 理 中 心 ， 各 台 站 的 antabfs 文件 整合 在 一 起 ， 生 成 包含 所 有 台 站 系统 温度 和 增 
益 信息 的 antab 文件 ， 供 数据 处 理 人 员 下 载 使 用 。 由 于 我 们 观测 时 的 部 分 台 站 没有 系统 温度 信息 ， 因 
此 在 AIPS 中 进行 幅度 校准 时 ， 需 要 手动 编辑 antab 文件 。 我 们 可 以 根据 每 个 天 线 的 标 称 系统 温度 计 
算 单位 流量 密度 ( Degree Per Flux Unit, DPFU) ， 编 辑 antab 文件 ， 然 后 再 通过 AIPS 任务 APCAL， 手 
动 进行 幅度 校准 。 
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Clean RR map. Array: BHKSSZ Clean RR map. Array: HIKKKKPSTTUY 
CTA 26 at 8.604 GHz 2018 Jul 20 1454-060 at 8.391 GHz 2016 Jul 28 
I I 


Relative declination / mas 
Relative declination / mas 


10 5 0 =5 -10 10 5 0 =) =10 
Right ascension / mas Right ascension / mas 
Map center: RA: 03 39 30.938, Dec: -01 46 35.804 (2000.0) Map center: RA: 14 56 41.393, Dec: -06 17 43.260 (2000.0) 
Map peak: 0.745 Jy / beam Map peak: 0.038 Jy / beam 
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Al 参考 源 图 像 。(a)HR 1099 HAAR CTA 26; (b)Del Lib 的 参考 源 1454-060。 图 中 左下 角 椭 圆 为 图 像 的 综合 波束 
Fig. 1 Images of the calibrators. (a) The calibrator CTA 26 of HR 1099; (b) the calibrator 1454-060 of Del Lib; 


the synthesized beam is shown in the bottom left corner of each image 


我 们 从 IVS 或 LBA 的 网 站 下 载 了 部 分 台 站 的 antab 文件 ， 对 于 缺少 信息 的 台 站 ， 我 们 手动 添加 
形 如 “081 08; 47: 52 600 600” 的 文本 到 antab 文件 中 ， 即 “观测 当天 在 该 年 中 的 天 数 - 当 天 时 间 
(UTC) -左旋 系统 温度 - 右 旋 系 统 温度 ”( 本 文中 仅 编辑 右 旋 数 据 ) 。 系 统 温度 7 可 以 通过 公式 SEFD- 
T../DPFU 计算 ,各 个 天 线 的 标 称 SEFD( 或 7 ) 可 以 通过 查询 IVS 和 LBA 官网 得 到 。 在 APS 中 导入 
antab 文件 后 ， 将 antab 文件 中 单一 时 刻 的 系统 温度 通过 APCAL 任务 应 用 到 完整 的 观测 时 段 。 将 AIPS 
校准 完成 后 的 参考 源 数 据 导入 DIFMAP, sos TI RE RSME, i51] gscale 命令 后 ， 根 据 结 果 ， 再 次 修 
改 刚才 步骤 中 的 antab 文件 ， 运行 gscale 后 结果 的 值 是 各 台 站 每 通道 名 义 SEFD( 根据 antab 文件 里 的 
DPFU 和 “FI7” 值 计算 得 到 ) 的 平方 根 ， 若 要 改正 antab 文件 中 的 “F7” 值 ， 需 要 对 gscale 结果 中 各 通 
道 值 分 别 求 平方 ， 取 中 位 数 后 更 新 到 “F7” 值 上 ， 从 而 可 以 将 gscale 的 值 修正 到 0.9~1.1， 以 此 校准 

台 站 的 幅度 信息 。 此 步骤 经 过 2~4 次 迭代 ， 可 以 将 各 个 台 站 gscale 值 修正 到 1 附近 ， 至 此 幅度 校 
bp Ars 述 幅 度 校准 过 程 中 各 人 台 站 系统 温度 信息 仅 为 大 致 估算 ， 所 以 不 能 精确 得 到 目标 源 的 流量 信 
息 ， 但 具有 一 定 的 参考 价值 。 

我 们 分 别 将 AIPS 校准 完毕 后 的 文件 导入 DIFMAP ， 对 两 颗 射 电 星 进行 成 图 ， 使 用 自然 权 得 到 两 

颗 射 电 星 的 图 像 。 


3 ”结果 分 析 及 讨论 


3. 1 图像 结果 分 析 
图 2(a) 为 射电 星 Del Lib AY IVS 搭载 式 相 位 参考 观测 的 图 像 结果 ， 由 于 射电 星 Del Lib 的 源 跟踪 时 
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间 仅 有 36 min AA, BNA 4 h, TUARUA F f8 BU UR EROR, BAS Se, 
UV 覆盖 东西 方向 较 差 ， 南 北方 向 较 好 ， 导 致 东西 方向 分 辩 率 较 差 ， 所 以 图 像 的 旁 汶 效应 很 明显 。 表 
3 列 出 了 Del Lib 的 单 历 元 坐标 和 误差 ， 表 中 亦 纬 方向 形式 误差 较 大 。 文 [8] 的 观测 显示 ， 该 射电 星 在 
8.3 GHz 的 宁静 阶段 流量 约 为 1~ 10 mJy， 爆 发 阶段 的 光斑 流量 达 20 mJy， 最 高 达到 34 mJy。 我 们 在 
8.4 GHz 频率 通过 粗略 antab 幅度 校准 测 得 的 流量 约 为 1 mJy， 图 像 光 束 尺寸 为 2. 38 x 0.51 mas, 


0 1 m Jy / beam 0 1 m Jy / beam 


00 35 12.952 
-08 31 08.352 


12.949 

3 g 

8 08.358 sj 

2 H 

a ja) 
12.946 

08.364 
15 00 58.2786 58.2782 58.2778 03 36 47.2528 47.2526 47.2524 
RA (J2000) RA (J2000) 
(a) (b) 


图 2 #22 AR, (a) A Del Lib, (b): HR 1099; 轮廓 开始 于 3 倍 的 图 像 均 方 根 ， 增 长 因子 为 2; 
图 中 左下 角 椭 圆 为 图 像 的 综合 波束 
Fig.2 Images of the radio stars. (a) Del Lib, (b) HR 1099; contours start at three times of the noise level of images 


and increased by factors of 2; the synthesized beam is shown in the bottom left corner of each image 


图 2(b) 为 射电 星 HR1099 的 图 像 结 果 ， 可 以 看 到 一 个 明显 的 左边 成 分 较 大 、 右 边 成 分 较 小 的 相 
接 结构 。 文 [22] 和 文 [23] 提出 可 以 用 核 晕 (core-halo) 结构 描述 。 文 [12] 指 出 ，HR 1099 的 VLBI 图 像 
结构 呈现 左 大 右 小 的 两 个 相 接 在 一 起 的 成 分 ， 呈 东西 方向 ， 角 距 约 为 1.7 mas， 他 们 在 形态 学 上 将 这 
个 “envelope” 称 为 “双核 坚 ” 或 “ 核 咖 ”结构 。 我 们 的 VLBI 图 像 结果 也 呈现 “ 核 党 ”结构 ， 呈 东西 方 
向 ， 角 距 约 为 2 mas， 与 前 人 结果 相符 。 前 人 的 观测 表明 该 系统 的 流量 非常 强 且 变化 很 大 。 该 射电 星 
在 宁静 阶段 的 流量 水 平 在 5~25 mJy(4. 86 GHz) 0 。 在 爆发 阶段 ，HR109%9 的 射电 流量 在 数 小 时 内 发 
生 显 著 和 快速 的 变化 ， 有 两 次 在 1 天 内 从 大 约 100 mJy 增加 到 1 Jy 775! RS CVn 双星 系统 的 射电 发 射 
通常 认为 是 来 自 延 展 磁 场 结构 中 的 弱 相 对 论 电 子 的 同步 辐射 。 经 antab 粗略 幅度 校准 后 ， 我 们 观测 
的 HR1099 在 8. 6 GHz 频率 的 流量 大 约 为 5. 4 mJy， 也 与 前 人 的 观测 相符 "站 ， 图 像 光 束 尺寸 为 0. 83 x 
0.71 mas, HR 1099 的 源 跟 踪 时 间 约 为 75 min， 总 的 时 间 跨 度 约 有 10h，UV 覆盖 较 好 ， 为 成 功 探 测 到 
“ 核 曙 ”结构 提供 了 有 利 条 件 。 

两 颗 射 电 星 的 亮 温度 可 由 公式 


Sot 1+ (1) 
ros 9 


Vas UU, 
HE, He, S, FEU EUER Dd, FAY mJy; zw 是 观测 频率 ， 单 位 为 GHz; 9, 和 9, 是 目标 
图 像 高 斯 椭圆 模型 的 半 功 率 全 宽 的 长 轴 和 短 轴 ， 单 位 为 mas; z 是 红 移 ， 因 为 目标 在 银河 系 内 ， 所 以 
此 处 为 0。 由 (1) 式 佑 算得 到 Del Lib 的 亮 温 度 约 为 2 x 10' K, HR 1099 的 亮 温度 约 为 5x 107 区， 两 者 
均 为 非 热 辐射 ”|。 


T, =1.22 x 10° 
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3.2 坐标 结果 分 析 及 视差 月 行 拟 合 
表 3 给 出 了 两 颗 射 电 星 单 历 元 坐标 和 图 像 信 噪 比 等 信息 ， 其 中 坐标 误差 是 由 JMFIT 给 出 的 形式 误 
差 。 表 3 最 后 一 列 cwiw 为 射电 星 的 整体 定位 精度 (位 置 误差 ) ， 可 由 经 验 公 式 


C. bosition = E J i d (2) 
T) J8ln2 5 
Ah. Hob, 0 是 波束 宽度 的 半 功 率 全 宽 ; oy 为 流量 的 均 方 根 ; $ 为 图 像 上 流量 峰值 ， 也 即 C 
Tu 1 ; 0 
= [ee JAM ——. Jl 经 过 化 简 JL ES S o 80.5 — 
信 噪 比 的 倒数 TTE 则 (2) 式 经 过 化 简 后 可 以 写 为 oi6, 守 0.5 SNR? 


表 3 同时 给 出 了 前 人 得 到 的 这 两 颗 射 电 星 单 历 元 坐标 和 精度 。 对 比 可 以 发 现 ， 我 们 的 精度 要 比 文 
[31-32] 的 高 一 个 量 级 ， 原 因 可 能 是 由 于 我 们 观测 的 基线 相 比 于 VLA + PT( Very Large Array + Pie Town) 
要 长 得 多 。 

表 3 目标 射电 星 的 坐标 等 信息 
Table 3 The coordinates of the radio stars and other information 


Star R. A. Dec. Oa Os Date Project code/ Flux Beam size — C position 


name (hms) (9: /mas /mas / year Data from j /m]y / mas / mas 


Del Lib 15 00 58.2781083 -08 31 08. 356379 0.07 0.26  2016.57203 | AOVOIO! 9.3 1.0 2.38x0.51 0.08 


15 00 58. 3456 -08 31 08. 219 13 18 2000.9426 D. Bobol 37.5 1.4 - = 
15 00 58. 3328 -08 31 08. 248 7 11 2003.6933 D. Boboltz? 
HR1099 03 36 47. 2525923 +0035 12. 948236 0.07 0.05  2018.54817  VSISC! 14.0 5.4 0.83x0.71 0.03 
03 36 47. 2869 +00 35 15. 772 5 5 2000.9426 D. Bobolz 204.1 13.5 - = 
03 36 47. 2457988 +0035 12. 688260 0.02 0.04 2020.01545 S. Lunz4 
注 :! 数 据 来 自 本 文 ， 目 标 射电 星 的 绝对 坐标 及 误差 由 AIPS 任务 JMFIT 得 到 ， 误 差 仅 为 形式 误差 ; 信 噪 比 和 流量 由 DIFMAP 得 到 ; JOR 


宽度 为 图 像 的 综合 波束 大 小 ; ?数据 来 自 文 [ 31] 的 VLA+PT 观测 ; ?数据 来 自 文 [32] 的 VLA+PT 观测 ;“ 数 据 来 自 文 [33] 的 VLBA 观测 。 


表 4 给 出 了 我 们 结合 前 人 数据 使 用 最 小 二 乘法 拟 合 得 到 的 自行 和 视差 ， 表 中 还 列 出 了 Gaia 
DR30 和 文 [31-32] 的 结果 作为 对 比 ， 后 者 由 于 数据 量 不 足 ， 仅 有 自行 结果 。 尽 管 文 [ 31-32] 的 数据 
为 多 历 元 观测 ， 但 并 没有 给 出 具体 每 一 个 历 元 的 数据 ， 因 此 我 们 只 能 使 用 文中 的 单一 历 元 数据 进行 视 
差 拟 合 。 


表 4 自行 和 视差 拟 合 结果 及 对 比 


Table 4 The fitting results of proper motions and parallaxes and comparison 


Star name TI/mas 11, / ( mas * year ! ) u,/ (mas + year! ) Data from 
Del Lib 10. 101 + 8. 006 -63. 671 + 0. 419 -4.922 + 1.021 AOVO010+D. Boboltz?' ^? 
一 -66. 16 € 1.51 -6. 19 x 1. 62 D. Boboliz!*”? 
9. 282 x 0. 473 -63. 051 + 0. 385 —6. 024 + 0. 326 Gaia DR3!! 
HR 1099 — 34.6110. 174 -31.677 € 0. 076 -161.628 0.040 ^ V515C- D. Boboltz ^" € S. Lunz 
— -31.87 +0.41 -161. 09 + 0. 52 D. Boboltz ^" 
33. 753 +0. 087 -32. 894 + 0. 131 -161. 772 € 0. 118 Gaia DR3!" 


MAA VL, ATA w= pcos6 Fl u, us 


我 们 使 用 3 个 点 (AOV010 加 上 D. Boboltz 的 2 个 点 ) 对 Del Lib 进行 视差 自行 的 拟 合 ， 结 果 显 示 ， 
我 们 拟 合 的 Del Lib 自行 在 * 方向 与 Gaia DR3 结果 和 精度 都 非常 吻合 ，7y 方向 结果 略 有 差异 ， 但 也 在 
误差 范围 内 ; Del Lib 的 视差 IL 误差 较 大 ， 可 能 是 由 于 文 [31-32] 的 数据 精度 较 低 导致 的 ; 我 们 的 Del 
Lib 拟 合 结果 的 精度 均 高 于 文 [31-32]。 
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我 们 使 用 3 个 点 (V515C 加 上 文 [31-32] 和 文 [33] 各 1 个 点 ) 对 Del Lib 进行 视差 自行 的 拟 合 。 结 
果 显 示 我 们 拟 合 的 HR 1099 自行 与 文 [31-32] 和 Gaia DR3'" 的 结果 都 相符 ， 且 精度 要 高 一 个 量 级 ; 
视差 TT 在 数值 上 与 Gaia DR3'' 相符 ， 但 误差 要 高 一 倍 。 


4 总 结 和 展望 


我 们 利用 测 地 模式 中 搭载 的 相位 参考 观测 ， 使 用 不 长 的 观测 时 间 得 到 了 目标 射电 星 的 VLBI 图 
像 、 流 量 和 单 历 元 坐标 ， 然 后 结合 前 人 的 数据 ， 使 用 最 小 二 乘法 拟 合 得 到 了 它们 的 自行 和 视差 ， 与 前 
人 的 结果 相符 ， 并 且 提 高 了 一 定 的 精度 。 由 于 测 地 观测 仅 有 右 旋 数据 ， 因 此 ， 如 果 为 了 人 研究 射电 星 的 
左右 旋 偏 振 等 其 他 性 质 ， 则 需要 正式 的 射电 星 VLBI 天 体 测量 观测 。 

本 文 的 工作 证 明了 利用 测 地 观测 中 搭载 式 相 位 参考 模式 观测 射电 星 的 可 行 性 ， 未 来 可 以 利用 这 种 
方式 ， 将 此 模式 作为 正式 的 射电 星 VLBI 天 体 测量 观测 的 先导 试验 ， 验 证 所 选 射 电 星 的 流量 强度 可 否 
被 观测 到 ; 也 可 以 在 大 量 的 测 地 观测 中 穿插 此 类 搭载 式 射 电 星 观测 ， 以 充分 利用 观测 时 间 ， 作 为 射电 
HE VLBI 天 体 测 量 的 补充 ， 丰 富 射 电 星 VLBI 观测 的 数据 样本 ， 为 光学 和 射电 参考 架 连 接 做 贡献 。 

Bb. 感谢 上 海天 文 台 VLBI 相关 处 理 中 心 和 澳大利亚 LBA 相关 处 理 中 心 在 数据 相关 处 理 方 面 的 
协助 。 
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Abstract; The GCRF ( Gaia Celestial Reference Frame) based on Gaia astrometric data has systematic 
differences from the ICRF ( International Celestial Reference Frame) based on VLBI ( Very Long Baseline 
Interferometry ). An external independent evaluation of the data from Gaia is required to link optical and radio 
reference frame. Celestial bodies with both optical and radio radiation are generally used to link the two 
frames, quasars are used at the dark end (magnitude Gz:15) , and radio stars are used at the bright end (GS 
13). It is very difficult to apply for the ordinary VLBI observations, so we try to focus on some geodesic 
observations that can use phase reference mode to observe radio stars in geodesic observations. In order to test 
the feasibility of this kind of observation, we use piggy-back phase referencing mode in geodetic observation of 
IVS (The International VLBI Service for Geodesy and Astrometry) and LBA (Large Baseline Array), and 
select Del Lib and HR 1099 as the targets. Our observations have obtained the VLBI images, flux densities and 
single-epoch coordinates of these two radio stars, and the coordinate accuracy is higher than that of the 
previous work; at the same time, we also combine the data of the previous work to fit parallax and proper 
motion of the two radio stars. Our results are consistent with the results of the previous VLBI work and Gaia's, 
and some of the accuracy is higher. Our work verify the feasibility of observing radio stars with piggy-back 
phase referencing mode in geodetic observation, and lay a certain foundation for the subsequent enrichment of 
VLBI astrometric samples of radio stars and linking of optical and radio reference frames. 
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